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Abstract 

The mean-eddy interaction in a flow is investigated in terms of the energy transfer between 
its mean and eddy processes. In the Reynolds decomposition framework, the classical energetics 
formalism does not have the transfer faithfully represented: Energy transferred from the mean 
process to the eddy process is not equal in amount to the energy transferred in the opposite 
direction. As a result, the classical "Reynolds stress extraction", TZ = — v'v' : Vv, cannot be used 
to measure the mean-eddy interaction, or turbulence production/laminarization in the context of 
turbulence research. This paper shows that a faithful formalism can be rigorously established 
within the same framework. The resulting transfer sums to zero everywhere over the decomposed 
subspaces, representing a mere redistribution of energy between the mean and eddy processes, 
without generating or destroying energy as a whole. For this reason, it has been referred to as 
perfect transfer in distinction from other energy transfers one may have encountered. The perfect 
transfer has a simple form in expression, with the introduction of an eddy flow coupled with the 
mean and eddy processes of the field under concern. The eddy or turbulence production is then 
totally determined by the divergence of this flow. This formalism has been validated with a well- 
known barotropic instability model, the Kuo model for the stability of the zonal atmospheric jet 
stream. We see a distribution of perfect transfer consistent with the instability scenario inferred 
based on Kuo's theorem, while the traditional Reynolds stress extraction TZ does not agree with 
the inference. The formalism has also been validated with a wake control problem. It is found that 
placement of a control to defeat the positive perfect transfer T (mean-to-eddy transfer) will yield the 
most efficient control in terms of energy saving, in comparison to the inefficient control placement 
based on TZ or eddy energy growth. This research is expected to be useful for the harnessing of 
turbulence, in that it may help to identify the best locations to place passive controls, or design 
the performance functional for active controls. 
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I. INTRODUCTION 



The mean-eddy interaction in fluid flows is an important problem in fluid mechanics. 
Related to it are hydrodynamic stability, turbulence production, laminar izat ion, atmospheric 
cyclogenesis, hurricane generation, ocean eddy shedding, to name but a few. Central to the 
problem is the transfer of energy between the mean and eddy processes as decomposed 
(cf. Fig. [T]). The purpose of this paper is to quantify this transfer within the traditional 
Reynolds decomposition framework, and use it to investigate a new strategy of fluid control. 
In a forthcoming paper, this formalism will be extended to a more generic framework for real- 
time problems (Liang et al., manuscript submitted to SIAM J. Multiscale Model. Simul.) 




FIG. 1: A schematic of the mean-eddy interaction, which is characterized by the energy transfer 
F between the mean and eddy processes. 



The classical formalism of energy transfer can be best illustrated with the Reynolds 
decomposed equations for the advection of a scalar field T = T -|- T' in an incompressible 
flow V, where the overbar stands for an ensemble mean, and the prime for the departure 
from the mean. In the absence of diffusion, T evolves as 

^ + V-(vT) = 0, (1) 

whose decomposed equations are 



0, (2a) 



v'T' - v'T') = 0. (2b) 

Multiplying fl2a|) by T, and fl2b|) by T', and taking the mean, one arrives at the evolutions 
of the mean energy and eddy energy (variance) (H ^ 

TV • (VT') (3a) 



v'T' ■ VT. (3b) 

The terms in divergence form are generally understood as the transports of the mean and 
eddy energies, and those on the right hand side as the respective energy transfers. The latter 



dT - 

— + V ■ (vT + v'T') 



or 

'dt 



+ V ■ (v'T + vT' 



aTV2 
aT'V2 



+ V • (vT72) 



+ V • (vT'72) 
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are usually used to explain the mean-eddy interaction. Particularly, when T is a velocity 
component, the right hand side of (13bp has been interpreted as the rate of energy extracted 
by Reynolds stress, or "Reynolds stress extraction" for short, against the mean field to fuel 
the eddy growth; in the context of turbulence research, it is also referred to as the "rate of 
turbulence production" . 

An observation of the two "transfer terms" on the right hand sides of (I3D is that they 
are not symmetric; in other words, they do not cancel out each other. In fact, they sum to 
V ■ (Tv'T'), which in general does not vanish. This is not what one expects, as physically 
a transfer process should be a mere redistribution of energy between the mean and eddy 
processes, without destroying or generating energy as a whole. These two quantities therefore 
are not real transfers, and cannot be used to measure the mean-eddy interaction. 

The reason for the asymmetry between the terms on the right hand side of ([3]) is that 
they are intertwined with transport processes; or alternatively, the divergence terms on the 
left hand side do not account for all the fluxes. Some people such as Pope[3| add an extra 
term in the flux term of fl5al) to maintain the balance, but it is not clear how that term 
should be chosen on physical grounds. Pedloskyj^ pointed out that a partial solution of 
the problem is to take averages for these terms over a substantially large domain. This 
way the transport contributions may be reduced and hence the transfer stands out. Liang 
and Robinson [5| argued that spatial averages should be avoided to retain the information 
of spatial intermittency in the energetics. They believed that a precise separation between 
the transport and transfer can be made to satisfy the above symmetric requirement. They 
even named the transfer thus obtained perfect transfer, in distinction from other transfers 
that may have been called in the literature. But in their paper a rigorous formalization was 
postponed to future work; how the separation can be achieved is still open. 

This study intends to give this problem a solution in the traditional framework. A 
complete answer to the issue of separation raised in [5(] , which is based on a new mathematical 
apparatus, the multiscale window transform developed by Liang and Anderson (manuscript 
submitted to SIAM J. Multiscale Model. SimuL), is deferred to the sequel to this paper. 
The following two sections are devoted to the establishment of a rigorous formalism for 
the transfer F. We first consider the case for a scalar field T ([11]), and then extend to 
momentum equations fillip . The formalism is validated with a well-studied instability model 
of an atmospheric jet stream (section [TVT) . and applied to harness the Karman vortex street 
behind a circular cylinder (jV]). A brief summary is presented in section IVTl 



II. FORMALISM WITH A PASSIVE SCALAR 
A. Reynolds decomposition 

The transfer is sought within the Reynolds decomposition framework. The key of the 
Reynolds decomposition is Reynolds averaging. It decomposes a field, for example a scalar 
field T, into a mean T plus a departure from the mean, T' . Simple as it is, Reynolds averaging 
actually introduces an important geometric structure which, as we will see shortly, helps to 
make the transfer problem easier. 

A Reynolds average may be understood either as an ensemble mean, or an expectation 
with respect to the measure of probability. Practically it may also be understood as an 
average in time or an average in some dimension of space. Its basic properties include: 
(1) T' = 0, (2) aT = aT, for a = const, and a corollary from the above two, (3) TT' = 
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T(T') = 0. To put all these understandings together, the decomposition can be recast in 
the framework of a Hilbert space Ti, with an inner product defined as, 

{f,9)=79, (4) 

for any / and g in the space H, that is, the ensemble, the probability space, or the space of 
functions over the time or spatial domain under consideration, if the overbar is, respectively, 
an ensemble mean, a probability expectation, or a time/spatial average. (It is interesting 
to note that the meaning of a Reynolds averaging is two-fold: one is the mean state recon- 
struction, another the summation or integration operator in forming the inner product.) 

A Reynolds decomposition thus splits Ti. into two subspaces, which contain the mean 
process and the eddy process. We will refer to these subspaces as windows, so we have a 
mean window and an eddy window. Distinguish them respectively with subscripts and 1. 
Correspondingly the decomposed components of a field T, T and T', will be alternatively 
written as [T]q and [T]-^ for convenience. Using these notations, the energy of T on window 

k as defined in ([3]) is [T]^); the property TT' = becomes [T]-^) = 0, implying 

the two windows are orthogonal. The concept of orthogonal windows puts the mean and 
eddy fields on the same footing, and will help to greatly simplify the derivation. 



B. Multiscale flux 

An important step toward the solution of the transfer problem is finding the fluxes, and 
hence the transports, on the two scale windows. The (vT^/2) and (vT'2/2) in fl3aj) and 
(I3b|) . though seemingly in flux forms, are not really the desiderata in a rigorous physical 
sense. One may see this through a simple argument of energy conservation, which requires 
that the mean and eddy fluxes sum to (vT2/2)-clearly these two quantities do not meet the 
requirement. 

On the other hand, the concept of multiscale flux can be naturally introduced within 
the formalized Reynolds decomposition framework. Given a flow v, the flux of an inner 
product (/, g) over H is (v/, g) = (/, yjj) (x is self-adjoint with respect to (-, ■)). (Note 
the flux is uniquely represented this way. The only other choice one might propose for the 
representation is (v, fg). This, however, does not make sense in physics, as a flow v is an 
operator, not a function of the same like as / and g in the functional space.) Let f = 9 = 
one obtains the flux of energy 

Q = ^(vT, T). (5) 

Geometrically, the right hand side of shows a projection of (vT) onto T. The flux on 
window k then should be a projection of vT onto [T]^: 

Q, = ^(vT, [T],) = i([vT]„ [T],). (6) 

In arriving at the last result we have used the fact that the two windows are orthogonal. 

The multiscale fluxes thus obtained are additive, i.e., Q = + Q^- In fact, by the 
orthogonality between the mean and eddy windows, we immediately have 

E^([vT]„ [T],) = ^(vT, T). 

fe=0 

This is the very conservation requirement mentioned above. 
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C. Perfect transfer 



Continue to examine the evolution of T in an incompressible flow v. In the language 
introduced in subsection III At Eqs. ( !2al) and ( l2b1) can be written in a unified form: 



dt 



V • [vT], = 0, 



(7) 



for windows /c = 0, 1. Here the decomposition is performed only in a statistical sense. That 
is to say, the mean is an ensemble mean or a probability expectation. But as we will see 
toward the end of this section, the formalism of energy transfer is essentially the same with 
respect to other methods of averaging. 

Application of ([T]^, ■) to ([7]) gives the energy evolution equation on window k: 



dt 



+ {[T]u. V-[vT],)=0. 



(8) 



The nonlinear term (second part on the l.h.s) involves two interwoven processes: transport 
and transfer. The former integrates to zero over a closed spatial domain; the latter sums to 
zero over A;, A; = {0, 1}. That is to say, ([8]) can be symbolically written as. 



dt 



-v-Q 



(9) 



where is the flux on window /c, and the transfer to window k from its complementary 
subspace. 

We already know the multiscale flux in ([6]). The transfer Tk is now easy to derive. 
Comparing ([9]) with ([H]), one obtains 

Tk-^-% = -{[T],, V-[vT],). (10) 

Substitution of ([6]) immediately gives 



r. = 2V-([vT]„ [T],)-([r],, v-[vr],). 



(11) 



It would be more clear to see the mean-eddy interaction if ffTTl) is rewritten in the tra- 
ditional overbar/prime notation. For the eddy process [k = 1), the transfer from the mean 
flow is 



Ti = ■ (^(vT)'T') - T'V ■ (vT)', 



(12) 



which reduces to 



Fi = F = - {TV • (v'T') - (v'T') ■ VT} 



(13) 
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In the derivation, the incompressibihty assumption V-v = 0, and hence V- v = 0, V- v' = 0, 
has been used. Likewise, 



1. 



and 



To = 2 V ■ ((vT)T) - TV • (vT) = -F, 



1 



Qi = -|vT'2 + rv'T' 
The mean-eddy energetics corresponding to ([3]) are, therefore, 

-—^ + V • -vT^ + -Tv'T' 
9t V 2- 2 - 



1, 



(fr''^ /2 

, ' + V • ( -vT'2 + -Tv'T' 
dt V 2- 2 - 



(14) 

(15) 
(16) 

(17a) 
(17b) 



with r as shown in f[T^ . 

Equations ([T^ and ([T^ imply an important property for the transfer derived above. 




(18) 



That is to say, the transfer thus obtained is a process of energy redistribution between the 
windows; there is no energy generated or destroyed as a whole, just as one may expect. To 
distinguish from other energy transfers one may have encountered in the literature, we will 
refer to F as perfect transfer, a term adopted from when confusion might arise. 

Note the distinct difference between F and the Reynolds stress extraction as appears in 



fl3bp . TZ = — v'T' ■ VT, which traditionally has been used to interpret the generation of 
eddy events, and has been interpreted as, in the turbulence research context, the turbulence 
production rate. In sections [IV] and El we will see that these two are in general differently 
distributed in space and time. 

The transfer F may be further simplified in expression. If T 7^ 0, (fT3l) may be alternatively 
written as 



F 



v'T' 



(19) 



Observe that the quantity in the parenthesis has the dimension of velocity. It represents 
a flow coupled with the mean and eddy processes of T. For convenience, introduce a "T- 
coupled eddy velocity" 



v'T' 



(20) 



then 



;t'v ■ 



(21) 



6 



Notice that fT^ is the mean energy of T and is hence always positive, so whether eddies are 
produced is totally determined by the divergence of the T-coupled eddy flow. 

The T-coupled eddy flow v^. is introduced for notational simplicity and for physical 
understanding. One should be aware that V is well defined, even though v-^ does not exist 
when T = 0. In that case, the original expression (fT3l) should be used. 



D. Other methods of averaging 

Practically the Reynolds averaging is often performed with respect to time if the process 
is stationary, or some dimension of space if the process in homogeneous in that dimension. 
If the averaging is in time, the above derivations also apply, except that the time derivatives 
in ([7]), ([8]), fll7ap . and fll7bl) are gone. The transfer V is still in the same form as f|T3|) . 

If the averaging is performed in a dimension of space, say, in a;, then the above derivation 
needs modification, as the averaging does not commute with But we have the following 

extra properties: |^ = 0, ^ = 0, for any field ip. These substituted into the continuity 
equation V ■ v = yield |^ = |^ + ^ = 0. With these identities, we repeat the procedures 
in the above subsection, and obtain: 

r = ^ {TV,. ■ {VT) - {yT) ■ V,.f } . (22) 

Here V yz = + k.|^ is the V operator with the x component removed. 
Notice that u'T' and T are independent of x, viz 



^^f^-(^l^)^0. (23) 



2 \ dx dx J 

So we may add the left hand side of (!23|) to (l22l) to get 

r = - {TV ■ (w^V) - (VT) ■ Vf } , (24) 

which is precisely the same as (1131) in expression form. 

In a brief summary, we have derived the mean-eddy energy transfer for a passive scalar 
in an incompressible flow, which is "perfect" in the sense that the energy extracted from 
the mean is equal in amount to the energy released by the mean. The perfect transfer is 
invariant in expression form with averaging schemes. 



III. FORMALISM WITH MOMENTUM EQUATIONS 

We need to deal with momentum equations for the mean-eddy kinetic energy transfer. 
Consider an incompressible ideal flow v. It is governed by 

^ + V ■ (vv) = -VP, (25a) 



_v_ — w. 



0, (25b) 



In forming the multiscale energy equations, the pressure term only contributes to the trans- 
port, i.e., the resulting pressure work is in a divergence form, thanks to the incompressibility 
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assumption. So only the nonlinear terms require some thought in deriving the perfect trans- 
fer. 

In this case, the formalism with the momentum equations is then essentially the same as 
that with the evolution of a passive scalar, except that we now have three "scalars" for the 
three velocity components u, v, and w. If m 7^ 0, -u 7^ 0, and w 7^ 0, for each component, 
there is a "coupled eddy velocity" defined by (120|) . so we have v^, v^, and v^, which are 

^4^, and respectively. According to the preceding section, each velocity corresponds to 
a transfer as expressed by fl2T|) . The total kinetic energy transfer is thence the sum of all 
the three transfers: 

r = 1 {u'V ■ v„ + v'V ■ V, + «;2V • v^} . (26) 

As noted in subsection III C\ the above formula is well defined even when the mean velocity 
vanishes. In that case, one just needs to expand it to obtain: 

r = ^ {V- (Vv) ■ V- (VV) : Vv} , (27) 

where the second term in the curly braces is the very Reynolds stress extraction: 

TZ = -{Vv) : Vv. (28) 

Like those with a scalar field, these formulas stay invariant in form, no matter what an 
averaging scheme is adopted. 



IV. VALIDATION WITH AN INSTABILITY MODEL 

In this section, the above formalism of perfect transfer F is validated with an idealized 
instability model. We will also see through this concrete example how F differs from the 
classical Reynolds stress extraction against the basic profile. 

Consider a well-studied barotropic instability model, the Kuo model, for the instability of 
the zonal atmospheric jet stream[g| 0]. Liang and Robinsonjsl have constructed a particular 
solution with a highly localized structure which is ideal for our purpose here. In the following 
we briefly present this solution, and then calculate the transfer (!27|) . 

Choose a coordinate frame with x pointing eastward, y northward. The governing equa- 
tions for the Kuo model are the 2D version of (!25|) . but with a Coriolis force term /k A v (/ 
constant) on the left hand side. The domain is periodic in x, and limited within latitudes 
y = ±L, where a slip boundary condition f = is applied. As rotation makes no contri- 
bution to the energy evolution, the formulas established in section (1IIII) equally apply here, 
i.e., the Kuo model can be used for the validation. 

Assume a basic velocity profile (cf. Fig. [2^) 

U{y) = ^^maxCOS^ (^f) ' ^ ^' ^^^^ 

The background potential vorticity q has a meridional gradient (cf. Fig. [2)d) 

Qy = -Uyy = ^ ^TJl^Ujai,^ COS — , (30) 
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(a) 



(b) 









qy<0, li-c r<0 
qy<0, u-c r>0 

qy>0, u-c r>0 

qy<0, u-c r>0 
qy<0, u-c r<0 



u/ u 



max 



qy/ (Umax/L) 



FIG. 2: Configuration of tlie Kuo model, (a) The basic flow profile u = u{y). (b) The background 
potential vorticity. Marked are the two reflection points on the profile curve. 



which changes sign at y = ±-|, meeting the necessary condition for instability by Rayleigh's 
theorem {ibid). Decompose the flow as 

{u,v) = {u{y),0) + {u',v'), (31) 

and substitute back to the governing equations. Kuo considered only the initial stage of 
instability when the perturbation field {u', v') is very small. So the resulting equations can 
be linearized. Assuming a solution of the form 

K,t;') = (u(l/),%))e^'=(--*), (32) 

one obtains an eigenvalue problem 

^+('^_iA£, = o, (33) 

with boundary conditions 

{} = 0, at ?/ = ±L. (34) 

The solution of (l33l) is not repeated here; the reader may refer to Kuo's original papers for 
details. 

Kuo showed that, in addition to the qy inflection requirement, the difference {u — Cr) 
{cr = Re{c} the mode phase velocity) must be positively correlated with qy over [—L, L] in 
order for the perturbation to destabilize the jet. In other words, for an instability to occur, 
it requires that 

(1) qy change sign through y G [—L,L] (Rayleigh' theorem); 

(2) (m — Cr) and (jy be positively correlated over [— L, L] (Kuo's theorem). 
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Hence the zero points of {u — Cr) and qy are critical. We will validate our transfer formal- 
ism through examining the instability structures near these critical points. We choose a 
particular unstable mode (and hence a particular c^) to fulfill the objective. 

As shown in jsj, the wavenumber ^ = |f gives such a mode; it lies within the unstable 
regime as computed by Kuo[7]. In fact, if substituting back into the eigenvalue problem, 
one obtains, using the shooting methodjOi], 

c = cr + ici = (0.4504 + 0.0476i)Mmax, (35) 

yielding a positive growth rate kci > 0. Solved in the mean time is the corresponding 
eigenvector v, which substituted in (152]) and the governing equations give a solution of all 
the fields. The resulting phase speed Cr = 0.4504Mmax and the gradient of the basic potential 
vorticity qy give four critical values of y: 

uiy) -Cr = ^y = ±0.53L, , . 

% = -Uyy = ^y = ±0.50L. ^ ^ 

The four critical latitudes, as marked in Fig. partitions the y dimension into five 
distinct regimes characterized by different values oi K = qy{u — c^). For most of e [— L, Z], 
K > but the positivity is interrupted by two narrow strips near y = ±L/2, where < 0. 
This scenario has profound implications by Kuo's theorem. Although Kuo's theorem is 
stated in a global form, it should hold locally within the correlation scale. In the present 
example, that means one of the necessary conditions for barotropic instability is not met 
around the strips and so there should be no instability occurring there. 

(a) PERFECT TRANSFER (b) REYNOLDS STRESS EXTRACTION 

1 I , , 1 , , 1 1 I , , 1 , , 1 




FIG. 3: The barotropic energy transfer (scaled by u'^g^^/L) for the Kuo's model: (a) the perfect 

here by (f27l) : (b) the Reynolds stress extraction TZ, which 



transfer F, which is 



- du'v' 

' dy 



is equal to 



dy 



here. The averaging is taken with respect to x. 



Instability means a transfer of energy from the background to the perturbation field, 
namely, a positive F. Using the particular solution obtained above, we compute the transfer 
from (l27j) . We adopt a zonal averaging, i.e., averaging in x, to fulfill the decomposition. 
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This is because, (1) u itself does not have x-dependence and hence can be understood as an 
x-average, and (2) the solution is homogeneous in x due to the cyclic boundary condition. 
The computation is straightforward. The result is plotted in Fig. [3^1. Sure enough, T is not 
positive around the two narrow strips; in fact, there is a strong negative transfer, i.e., upscale 
or inverse transfer from the eddy window to the background. Moreover, the negative transfer 
is limited within two narrow regimes, just as one may expect by Kuo's theorem. In contrast, 
a different scenario is seen on the profile of the conventional Reynolds stress extraction TZ, 
which we plot in Fig. [3]d. TZ is nonnegative throughout [—L, L]; particularly, it is maximally 
positive over the narrow strip regimes, countering our foregoing intuitive argument. Through 
this example, our perfect transfer F results in a scenario agreeing well with the analytical 
result of the Kuo model, while the conventional Reynolds stress extraction TZ does not. 



V. APPLICATION TO THE SUPPRESSION OF EDDY SHEDDING 

A practical application of the above research is turbulence control. Turbulence control is 
a technique to manipulate turbulence growth to achieve the goal of drag reduction (cf. [13] 



and particularly the celebrated paper by Kim[ll[, and the references therein). How the 
current research may come to help is to provide a better object, i.e., the perfect transfer, to 
manipulate, in place of the growth of turbulence energy or eddy energy. 

The proposal is out of the concern of how to maximally take advantage of the processes 
of self-laminarization or relaminarization that may occur within a turbulent flow. In the 
interest of energy saving, suppression of the positive transfer F [cf. fl27|) ] is preferred to 
suppression of the eddy energy growth to inhibit the production of turbulence. To see 
why, observe that eddy energy increase does not necessarily occur in accordance with a 
positive transfer, and hence a place where turbulence grows does not necessarily correspond 
to turbulence production. Actually, the correspondence is an exception rather than a rule. 
(Later in this section we will see an example.) The energy needed to fuel the growth could 
be transported from the neighborhood, rather than released in situ. One possibility is, 
while disturbances rapidly grow, a process of laminarization could be undergoing at the 
very position. As shown in the two-point system in Fig. HJ while disturbances grow at both 
A and B (both K'^'^^ and K'j^'^^ increase), the eddy energy is produced at A only. At B, 
not only there is no eddy energy production, but the transfer is from the eddy window to 
the mean window. That is to say, the system is undergoing a laminarization at B, even 
though the eddy energy K'^'^^ grows, because of a surplus of the influx of eddy energy over 
the inverse transfer. Control of the eddy energy growth at both A and B indeed helps to 
suppress the onset of turbulence, but it is not optimal in terms of energy saving. Suppression 
of K'^^^ defeats the intrinsic trend of laminarization in the mean time, and therefore reduces 
the control performance. To take advantage of this laminarization, the control should be 
applied at position A, i.e., the source region, only, and the optimal objective functional 
should be designed with respect to F^, rather than K'^'^^ + K^^^ . 

In this spirit, we demonstrate the application by showing how one may efficiently suppress 
the vortex shedding behind a circular cylinder. Vortex suppression is important in that it can 
result in significant drag reduction and hence energy saving; it may also be used to reduce 
noise. Presented in the following are just some diagnostic results with a saturated wake 
to which the afore-established formalism is applicable; the same example will be studied in 
more detail in the sequel to this paper with nonstationarity considered. We will deal with a 
laminar case only, but the idea equally applies to turbulent wakes. 
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FIG. 4: Schematic of the eddy energy transport and transfer for a two-point turbulent system. 
K stands for kinetic energy, and F and V • Q for transfer and transport, respectively. An arrow 
indicates the direction of an energy flow, with its thickness standing for strength. In this case, 
transfer is toward the mean at position B, but K^^^ still grows because of the transport (advection) 
of K'f'^y fr om position A. For optimal results, control should be placed at A only. 



There are many sophisticated techniques to suppress the shedding of vortices in a wake 
[cf. jHI and 13 1 and the references therein]. Surface-based suction is one of them. To our 
knowledge, the research along this line thus far, however, has been focused on the technique 
per se. No report has been found on the issue of where to place the suction to optimize the 
performance. In the following, we will show that our formalism of mean-eddy interaction 
and perfect transfer can give this question an answer. 

Consider a planar flow passing around a circular cylinder. The governing equations are 
the same as those of Eqs. fl2Sl) . except that dissipation is included. The computational 
domain is plotted in Fig. [5l with x and y nondimensionalized by the cylinder diameter d. A 
uniform inflow ([/, 0) is specified at a; = —2.5, and on the right open boundary (x = 37.5) 
a radiative condition [13] is applied, hi y = ±4 are two solid boundaries, where nonslip 
conditions are imposed. 

We examine a flow with Reynolds number Re = ^ = 200. The spacing choice of Ax and 
Ay is found not a stringent constraint. By experiments a mesh with Ax = Ay = 0.05 and a 
mesh with Ax = Ay = 0.025 produce little difference in the final result for our problem. We 
thus choose Ax = Ay = 0.05 for economic reason. The governing equations are integrated 
forward on a stag gered grid[l^ using a semi- implicit (implicit for pressure) finite difference 
scheme (e.g., [l6|) of the second order in both time and space. 

The model is run until a statistical equilibrium is reached. After that, it is integrated fur- 
ther for 100 time units and the outputs are used to calculate the transfer F. The stationarity 
in time makes it a natural choice to perform time averaging in computing the transfer fl27j) . 
The computation is straightforward. We plot the result in Fig. [6l Note the time interval is 
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FIG. 5: Model configuration. The coordinates x and y are scaled by the cylinder diameter d. 



large enough that one virtually sees no difference in the computed result if it is enlarged. 

Perfect energy transfer 




-4e-1 



12 3 4 

FIG. 6: Perfect transfer F in the wake behind a cylinder (units in U^/d). 



By ( |27I) . a positive F means eddy energy generation or turbulence production in turbu- 
lence research, while a negative F indicates a transfer in the opposite direction. In Fig. [6l 
two triangular lobes of strong positive F sit on either side of the axis y = 0, with a weak 
negative center lying in the near wake. That is to say, eddy energy is generated within 
the two lobes, while in between is a laminarization process. By the forgoing arguments, an 
efficient control strategy should be the one inhibiting the positive F in these two lobes. Since 
we are considering only the technique of surface-based suction, the F distribution suggests 
that application of suctions near the two lobes should be effective. In doing this, one can 
simultaneously take advantage of the laminarization process occurring in the near wake. 
Indeed, our control experiments show that the areas between 50 to 80 degrees and between 
—50 to —80 degrees from the x-axis are the effective suction locations to suppress the vortex 
street. The effectiveness, according to [l3|, may be measured by a suction rate Cg = 
where m is the mass flow rate, U the free stream velocity, and d the cylinder diameter. By 
experiments, the most effective control is that placed at ±70°, following the same orienta- 
tion of the two positive F lobes. Only a rate of Cg = 0.18 can have the vortices completely 
suppressed (e.g.. Fig. [7]). In contrast, controls in areas below 50 degrees and above —50 
degrees are counterproductive, as the near wake laminarization process is defeated. 
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t = 50 



t = 78 




FIG. 7: Snapshots of the vorticity in the optimal control experiment. The control is applied at 
t = 50 and so forth. 



It is of interest to see how other diagnostic fields, such as the eddy energy + f 
and TZ, are distributed. Shown in Figs. [Ht and b are these fields. Clearly, they both attain 
their maxima along the axis y = 0, a scenario completely different from that of F in Fig. El 
If the control is based on these fields, a suction should be placed at {x,y) = (|,0), i.e., 0° 
from the axis. But, as noted above, the control experiments show this does not result in an 
effective vortex suppression. 



(a) Eddy energy (Cl=0.05) (b) Conventional transfer 




1 2 3 4 1 2 3 4 



FIG. 8: (a) Perturbation energy (in U'^); (b) 7^ (in U^/d). 

The success of the control experiment serves to validate our formalism of the mean-eddy 
interaction and the equations of F, f|T3|) and fl271) . In the mean time, a variety of fluid control 
problems, both active and passive, may benefit from this formalism. 

VI. CONCLUSIONS AND DISCUSSION 

In the Reynolds decomposition framework, the mean-eddy interaction has been rigor- 
ously formulated in terms of energy transfer, which can be singled out from the intertwined 
nonlinear processes by eliminating the transport effect. The resulting transfer sums to zero 
over the two decomposed subapces, or windows as called in the text. In other words, the 
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transfer represents a redistribution process between the mean and eddy windows, without 
generating or destroying energy as a whole. Because of this property, it is sometimes referred 
to as perfect transfer in distinction from other transfers that one may have encountered in 
the hterature. 

The perfect transfer from the mean process to the eddy process can be exphcitly written 
out. In the case of a scalar T advected by an incompressible flow v, traditionally there is a 
quantity 

n = -{Vv) ■ vf , 

which, when T is a velocity component, has been explained as the rate of energy extracted 
by Reynolds stress against the basic profile. We showed that this is not the eddy energy 
transferred from the mean to the eddy windows. The real transfer should be 

r = i [f V ■ (yr) + n] , 

which may also be written as 

r = V v-v^, 

if r 7^ 0, in terms of a "T-coupled eddy flow" 

TV 

-T rp 

Since is the eddy energy and is hence always positive, the eddy generation is therefore 
completely controlled by the divergence of this flow. This simple formalism can be easily 
generalized to those with momentum equations. In that case, the perfect transfer is a 
redistribution of kinetic energy between the mean and the eddy windows. The resulting 
transfer is referred to fl27j) . For all the averaging schemes, it has the same form. 

The formalism has been validated with a well-known barotropic instability model, the 
Kuo model for the stability of a zonal atmospheric jet stream. Instability implies energy 
transfer from the background to perturbation, or mean to eddy in this context. We have 
seen a scenario of perfect transfer consistent with that inferred based on Kuo's theorem, 
while the traditional Reynolds stress extraction does not agree with the inference. 

An intuitive argument regarding the perfect transfer F is that the distribution of F is 
generally not in accordance to that of eddy energy or eddy energy growth, due to the 
presence of transport processes. This has been testified in the wake control experiment. 
In the context of turbulence, that is to say, the rapid growth of turbulent energy does not 
necessarily correspond to turbulence production. It is not uncommon that, at a location 
where perturbation is growing, the underlying process could be a transfer in the inverse 
direction, i.e., a laminarization. This argument has profound implication in real applications. 
Turbulence control is such an example. It suggests the optimal location to place a control 
be that of positive F, rather than that of turbulence growth, in order to take the advantage 
of the self-laminarization within a turbulent flow. This conjecture has been testified in an 
exercise of vortex shedding suppression with a cylinder wake. By computation there are 
two lobes of strong positive F attaching to the cyhnder on either side. We tried a surface- 
based suction on many places of the cylinder, but the most effective places are those where 
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the transfer processes within these two lobes are easiest to defeat. Other places are not as 
effective as these two, in terms of energy saving. 

The success of the wake suppression experiment implies the physically robust quantity F 
may be useful for a variety of fluid control problems. Speciflcally, it may come to help in 
selecting the location(s) to place a passive control, or designing the performance functional 
for an active control. The above experiment is an example for the former; for the latter, 
we should be able to design some transfer-oriented functional for the optimization. As we 
argued before, this should be advantageous over those based on turbulence growth in light 
of energy saving. 

It should be pointed out that, in realistic flows, the signals are generally not stationary, 
nor homogeneous, and as a result, the Reynolds averaging cannot be replaced with an 
averaging over time or a spatial dimension. In such cases, the mean and eddy flelds are not 
as simple as thus reconstructed; the mean itself can be time varying. Besides, interactions 
may not be limited just between two windows. A common process, mean-eddy-turbulence 
interaction, for example, requires three distinct windows for a faithful representation. All 
these difficulties will be overcome, and a new real problem-oriented formalism will be realized 
in a forthcoming paper after the introduction of a new analysis apparatus, multiscale window 
transform, to replace the Reynolds averaging technique for a realistic mean-eddy-turbulence 
decomposition. 
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